The Decay Properties of the Finite Temperature Density Matrix in Metals 
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Using ordinary Fourier analysis, the asymptotic decay be- 
havior of the density matrix F(r,r') is derived for the case 
of a metal at a finite electronic temperature. An oscillatory 
behavior which is damped exponentially with increasing dis- 
tance between r and r' is found. The decay rate is not only 
determined by the electronic temperature, but also by the 
Fermi energy. The theoretical predictions are confirmed by 
numerical simulations. 
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for the band-structure energy e n (k) and for the chemical 
potential ji. The integral extends over the Brillouin zone 
(BZ). 

The cell periodic part of the Bloch functions it n ,k(r) 
can be pulled out of the integral using the medium value 
theorem. Denoting the medium value by fco(r,r') one 
obtains 
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The decay behavior of the density matrix is a funda- 
mental property of quantum mechanical systems since it 
determines the degree of localitytl of all relevant quanti- 
ties. Experimentally it is well known, that most quanti- 
ties such as bonding properties and energies or magnetic 
moments depend only on an environment comprising a 
few nearest neighbor shells. Locality is certainly much 
stronger in insulators, but is also applies in metals. Re- 
cently there has also been renewed interest in these decay 
properties, since they determine whether and for what 
system size so called 0(N) methodsEl are faster than tra- 
ditional electronic structure algorithms. The decay be- 
havior of the density matrix in insulators has been the 
subject of several studies. It was shown that the decay 
constant is related to the location of branchunoints of the 
band structure in the complex energy planeou. Since this 
is a quantity which is not accessible in ordinary electronic 
structure calculations, it would of course be highly desir- 
able to relate the decay to features of the real part of the 
band structure. For the case of periodic solids arguments 
can therefore be found proposing a relation between the 
size of the gap and the decay behaviorETLl. Based on an 
analysis of the expansion of the density matrix in terms 
of Chebychev polynomials, this was recently confirmed!!! 
for arbitrary systems. 

In this paper the decay properties of the finite temper- 
ature density matrix in metals will be studied. Contrary 
to a widespread belief that the decay constant is given 
by y/hsT where ks is Boltzmann's constant and T the 
electronic temperature, we find that the decay is as well 
determined by the Fermi level. n 

In the finite temperature formalism of Mermincl the 
finite temperature density matrix F is given by 
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where ^ rl .k(r) = !in,k(r)e lk ' r ' are the Bloch functions 
associated with the wave vector k and n is the band in- 
dex. The function /( e "l j-" fi ) is the Fermi distribution 
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The part we have taken out is periodic in both k and r 
and bounded from above since the wavefunctions ti ra ,k(r) 
are normalizable. If n refers to the lowest band the ex- 
pression can also never vanish since u n ,k(jc) has no nodes. 
Hence the expression is bounded from below as well and 
it can therefore not influence the asymptotic behavior. In 
the case of higher bands the expression can in principle 
vanish and might lead thus a faster asymptotic decay. If 
the u n ^(r) has for instance an p z like character, the cell 
periodic part can identically vanish for certain directions 
of r and for certain values of r'. The expression can how- 
ever certainly not vanish for all values of r and r'. Thus 
the conclusion that the cell periodic part does not change 
the overall asymptotic behavior remains valid. 

In addition we will assume that the band structure can 
be approximated by the free electron band structure. The 
decay behavior of the metallic system is thus the same 
as the one of jcllium 
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Using elementary calculus this three-dimensional integral 
can be transformed into an one-dimensional one. Using 
integration by parts the Fermi distribution can be re- 
placed by its derivative /' and one obtains 



F(s) = 
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where s = |r — r'|. 

To proceed further with an analytic evaluation of the 
resulting expression one has to approximate /' by a re- 
placement function. The basic idea is that as the tem- 
perature tends to zero the part coming from /' becomes 
a delta function. So in the following steps we will replace 
/' by another representation of the delta function and we 
will use the rules which are valid for delta function for 
our sharply peaked finite temperature function as well 
even though it becomes a delta function only in the limit 
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of zero temperature. A good choice for the replacement 
function is the following 



/'(*) 
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sinh{bx) 
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where a = \ tan (|) and b = n tan (-|) . The coefficients 
a and b were chosen such that the replacement function 
has the same local and global behavior as the original 
function. The evaluation of the remaining integral then 
gives 
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cos(kps) 



(6) 



where c = 1 + \/2 and where kp is the Fermi wave- vector 
related to the chemical potential by hk F = u. 

So we see two prominent features in the behavior of 
the density matrix. It has oscillations whose wave- vector 
is given by the Fermi wave- vector kp and it is decay- 
ing exponentially, since for large values of r the expres- 
sion (l + \f2cosh (-c^V)) tends to exp(-c^f^r). 

Despite the approximative character of Equation ^ one 
obtains very good agreement in the jcllium case between 
the exact density matrix obtained numerically and its ap- 
proximation (Equation |6|) over many orders of magnitude 
as shown in Figure [l| 
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FIG. 1. Comparison of the exact decay behavior of the 
density matrix and its approximation of Equation |^. Three 
different thermal energies ksT of 1/4, 1/8 and 1/16 are plot- 
ted for Hf = 1. 



This confirms the well known fact that the amplitude 
decays algebraically in the zero temperature case. More 
specifically a decay proportional to r~ 2 is found. 

Since the derivation was based on the assumption that 
the band structure is free electron like, the results may 
not be applicable in the case of metals which a a very 
complicated Fermi surface arising from a band-structure 
which is completely different from the free electron band 
structure. 

I thank O. Gunnarson for his valuable comments, 
which considerably improved the paper. 
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Let us note that all the approximations used in the 
derivation become exact as the temperature tends to zero 
and consequently we obtain the correct zero temperature 
result. 

_, , Ankp ( . , . sin(skp)\ 

F( S ) = -—(cos( S kp) (7) 
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